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?3 ! Abstract 



Por the Edwards- Anderson Ising spin-glass model in three and four dimensions (Sd 
and Ad) we have performed high statistics Monte Carlo calculations of those free-energy 
O . barriers -Fg which are visible in the probability density Pj{q) of the Parisi overlap 

J-^ I parameter q. The calculations rely on the recently introduced multi-overlap algorithm. 

In both dimensions, within the limits of lattice sizes investigated, these barriers are 
K> , found to be non-self-averaging and the same is true for the autocorrelation times of 

^ \ our algorithm. Purther, we present evidence that barriers hidden in q dominate the 

canonical autocorrelation times. 



PACS. 75.10.Nr Spin-glass and other random models, 75.40.Mg Numerical simulation studies, 
75.50.Lk Spin glasses and other random magnets. 
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1 Introduction 

Spin glasses (for reviews see references |I|, 0, 0, ^) constitute an important class of materials 
whose low-temperature state is a frozen disordered one. In order to produce such a state, 
there must be randomness and frustration among the different interactions between the 
spins (magnetic moments). Frustration means that no single spin configuration is favored 
by all interactions. In real materials such competing interactions are for instance created 
by magnetic impurity moments. The study of spin glasses developed essentially since the 
middle of the 1970's and is based on three approaches: experiment, theory and computer 
simulation. 

Experimentally it is not hard to find spin glasses . One kind of widely studied systems 
consists of dilute solutions of transition metal magnetic impurities in noble hosts. The 
impurity moments produce a magnetic polarization of the host metal conduction electrons 
which is positive at some distances and negative at others. Because of the random placements 
of the impurities they have random, competing interactions with one another. Spin glass 
states have also been found in magnetic insulators and amorphous alloys. Properties 
analogous to those of spin glasses, with the electric dipole moment playing the role of 
the magnetic one, have been seen in ferroelectric-antiferroelectric mixtures. The universal 
behavior of the observed phenomena is a major reason for the interest in these systems. 

A freezing temperature Tc may be defined by a cusp in the ac susceptibility and has, for 
instance, been studied for Cu-0.9% Mn 0. Below this transition temperature characteristic 
non-equilibrium phenomena are observed. A typical experiment is the measurement of the 
remanent magnetization, see for a study of (Feo.i5Nio.85)75Pi6B6Al3. A spin-glass sample 
is rapidly cooled in a magnetic field to a temperature below the transition temperature 
and the observation is that the decay of the magnetization depends on the waiting time 
after which the field is switched off. This phenomenon is called aging and has also been 
found in other disordered or amorphous systems such as structural glasses, polymers, high- 
temperature superconductors, and charge- density wave systems. These large characteristic 
time scales suggest the presence of many equilibrium or metastable configurations with a 
distribution of free-energy barriers separating them. 

For free-energy barriers in spin glasses a major complication arises from the fact that 
there is no parametrization of the relevant configurations by a conventional thermodynamic 
variable. In his work on the mean-field theory of spin glasses Parisi generalized the 



concept of an order parameter. In later language []T], |^, ^, Q this is expressed as follows: A 
spin-glass realization is defined by a set of frozen, disordered exchange coupling constants 
J = {,Jik\ and for each realization the Parisi overlap parameter is defined by 

^^ 4 = 1 

where the sum goes over the total number A^ of spins of the system and the spin superscripts 
label two (real) replica of the same realization. For given J the probability density of q is 
denoted by Pj{q) and its cumulative distribution function is xj{q) = /^^ dq'Pj{q'). Average 
over the disorder defines the functions 

P{q) = [Pjiq)].. = ^Y.PM and x{q) = [xj{q\^ = J_^a;^(g) , 

where #J is the number of realizations considered. In the infinite volume limit below the 
freezing temperature an increasing continuous part of x{q) characterizes mean-field behavior 
of spin glasses, whereas in ferromagnets as well as in the droplet picture [§] of spin glasses 
x{q) is a step function. 

Analytical calculations in mean-field theory show that violations of the fluctuation- 
dissipation theorem in non-equilibrium dynamics determine the static function x{q) and vice 
versa 0, see |jTO| for a review. Numerical calculations in 3(i and Ad Ising spin glasses [0, |12| 
support that this relationship holds also in finite dimensions. Of course, the entire Pj{q) set 
contains more information than its mean P{q) (equivalently x{q)). In this paper we study 
the distribution of the minima in q of the Pj{q) probability densities. For given J' the non- 
trivial {i.e. away from q = ±1) minima are related to free-energy barriers of the disordered 
system J'. The other way round, it is presumably model dependent (and worthwhile to 
investigate) to what extent free-energy barriers of the system JT" are reflected in the minima 
of the Pj{q) probability density. 

Conventional, canonical Monte Carlo (MC) simulations do not allow for an efficient inves- 
tigation of the Pj{q) minima, because the likelihood to generate corresponding configurations 
in the Gibbs canonical ensemble is small. This problem is overcome by the multi-overlap MC 



algorithm |TE[ which samples with an uniform distribution in g. It belongs to the class of 



multicanonical and related algorithms [|T4[ |T5[, which allow to focus on rare configurations of 



the Gibbs ensemble. For instance, at first-order phase transitions in 3d, configurations with 
interfaces are suppressed according to exp(— aAmin), where a is the interface tension and 



Ajnin is the mininial area of the interface. For temperature driven transitions configurations 
with interfaces are found for E in the energy range Ei < E < E2 where E2 = Ei + AE and 
AE is the latent heat of the transition. To generate such configurations with a good statistics 
it is sufficient to sample with a weight factor w{E) ~ l/n{E), where n{E) is the spectral 
density. Similarly, interfaces for magnetic field driven first-order phase transitions can be 
generated by sampling with an appropriate weight function w{M) of the magnetization M 
of the sample. 

Once Pj{q) is determined, we define the associated free-energy barrier F| through the 
autocorrelation time of a Id Markov process which has the canonical Pj{q) distribution as 
equilibrium state. The barrier autocorrelation time r^ is then defined through the second 
largest eigenvalue of the transition matrix of this Markov process and the free-energy barrier 
is ln(r|). 

In previous literature [|1^, 0, |T8|, |I9|, ^, |21|] investigations of spin-glass barriers relied 



on various numerical and analytical methods, which are distinct from ours. The results of 
I6| , p!7| , p!8| , |T9|, ^ may be summarized as support of a scaling law F'j^'^ ~ A^^' ^ for canonical 



free-energy barriers in the mean-field limit below the freezing temperature. 

In the next section we describe our methods and give an overview of our MC statis- 
tics. Section ^ presents and interprets our numerical results for free-energy barrier in q. 
Conclusions and an outlook are given in the final section ^. 

2 Overview of Methods and Data 

The energy of the Edwards- Anderson Ising (EAI) [^ spin-glass model is given by 

E = -^Jik SiSk , (2) 

(ik) 

where the sum is over nearest-neighbor pairs of a (hyper) cubic lattice. The spins Sj as well 
as the coupling constants Jik take on the values ±1, with equal probabilities, i.e. the sum 
N'^ E(ifc) Jik is of order 1/y/N. 

In our calculations we combine the two copies (replica) of the same realization and 
simulate with a weight function 

w{q) = exp[-p{E^ + E^) + S{q)] . (3) 

Here /3 = Jo/ksT is the inverse temperature in natural units, E^ and E"^ are the energies 
of the respective replicas and S{q) has the meaning of the microcanonical entropy of the 



Parisi order parameter (|T]). The multi-overlap algorithm weights spin configurations with 
an overlap parameter q in such a way that a broad histogram in g, eventually covering the 
entire accessible range — 1 < 5' < 1, is obtained. This allows then for accurate calculations 
of the empirical probability density Pj{q) of the Parisi order parameter for realization J . 



Although an explicit order parameter does not exist, our simulation method ||T3[ is in this 



way similar to the multimagnetical [0, which for ferromagnetic systems is a very efficient 
way to sample configurations with interfaces. 

Our EAI simulations are performed on N = L'^ {d = 3,4) lattices at /5 = 1 {3d) 
and (3 = 0.6 {Ad). Both values correspond to temperatures T = 1/(3 below the freezing 
temperature of the respective model {(3^ = 0.90 ± 0.03 {3d) ^, (3^ = 0.485 ± 0.005 {Ad) 
2^). Table |1] summarizes the statistics we have assembled as well as the performance of our 
code. MC updates are given in units of sweeps. Our Jj^ realizations were drawn using the 
pseudo random number generators RANMAR [^] and RANLUX |2^ (luxury level 4). In 



the simulations themselves we always employed the RANMAR generator due to CPU time 
considerations. 

For each realization J' the simulation consisted of three steps: 

1. Construction of the weight function (Rf). Here we employed an improved variant of 



the accumulative stochastic iteration scheme discussed in Ref. |]2^, algorithmic details 



will be published elsewhere p8[. The iteration was stopped after at least 4 tunneling 
events 

(g = 0) ^ (g = ±1) and back (4) 

occurred. Our precise request was in 3d 10 tunneling events for L = 4, 6, and 8, and 
20 events for L = 12, but for a few cases with only 4 events requested. In Ad it was 
10 for L = 4, 20 for L = 6 and 20 to 30 for L = 8. In few cases, the system was 
tunneling so rarely between g = ±1 that we decided to abort the run and restart with 
a different random number seed, which in most cases led (eventually after multiple 
tries) to improved tunneling performance. After the weight function is constructed 
and kept fixed, the average number of sweeps it takes to create a tunneling event (^) 
defines the autocorrelation time of the multi-overlap algorithm which in the following 
is denoted by 



Of course, r™"^ depends on the realization J at hand, and on the parameters used in 
phase 1: random number seed, number of tunnehng events requested, etc.. 

2. Equihbration run. This run of ra x 65536 sweeps was done to equihbrate the system 
for given fixed weight factors (n = 1, 4, 16, 32 for 3(i L = 4, 6, 8, 12 and n = 2, 8, 16 for 
4(i L = 4, 6, 8, respectively). 

3. Production run. Each production run of data taking was concluded after at least 
20 tunneling events as defined in equation (^) were recorded. To allow for stan- 
dard reweighting in temperature we stored besides histograms of the Parisi overlap 
parameter also a time series of measurements for the order parameter, energies and 
magnetizations of the two replica. The number of sweeps between two successive 
points in a time series is adjusted in such a way that each time series is made of 65536 



measurements. This is done by an adaptive data compression routine [^. Together 
with the condition on the minimal number of tunneling events this ensures that the 
number of sweeps between two successive points in a time series is approximately 
proportional to r'^^^. Some re-weighting results were reported in references |]13| 
and 1^, publication of others is intended [^ . 



With each realization J we associate the free-energy barrier F% of the \d Metropolis- 



Markov chain [^ which has the canonical Pj{q) probability density as its equilibrium 
distribution. The transition probabilities Tj ., are given by 
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(6) 



where Wij {i 7^ j) is a probability a la Metropolis to jump from state q = Qj to q = q, 

(g, = ViV, ^e[-iV,-iV + 2,... + iV]), 



jKqj, 



(7) 



T fulfills the detailed balance condition (with Pj) and as a consequence it has only real 
eigenvalues. The largest eigenvalue (equal to one) is non-degenerate, and the second largest 
eigenvalue Ai determines the autocorrelation time of the chain, 
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Figure 1: Canonical Pj{q) and (flat) multi-overlap P™"'^(g) probability densities for our 
realization with the largest free-energy barrier in 3d {L = 8) and 4d (L = 8). 



and we define the associated free-energy barrier for realization J' as 



n 



In(r^) 



(9) 



For the simple double-peak situation of first-order phase transitions the autocorrelation time 
r^ is proportional to the ratio p^ax^pmm -^j-^gj^g 
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-^7 



PAqn 



max[Pj(g)] and Pf 



mm 

0<(jr<|qmax 



[pjm 



I.e. F| = ln[P^^^/PJ'^] + const. This leads in U to r| ~ e'"^™ and F% ~ aA^^^ + const., 
where Amin is the minimal area between the coexisting phase regions and a is the interfacial 
tension. Equation (^ is the appropriate generalization to a situation involving multiple 
minima and maxima. The autocorrelation time r^ has to be regarded as a lower limit to the 
canonical autocorrelation time t'^^^ for the Markov process where the spin variables are the 
dynamical degrees of freedom. The definition (^ takes only barriers in q into account but 
not other barriers which may well exist in the multi-dimensional configuration space. 

The matrix T in (^ is tri-diagonal and sign symmetric. This special form allows for 
easy calculation of all its eigenvalues ^^. The realizations with the largest thus obtained 
free-energy barriers in 3d and 4(i are depicted in figure 0. Both do not show a complicated 
landscape, but a plain two-peak structure. Besides the canonical Pj{q) probability densities 
the essentially flat probability densities P™"^(g) of the multi-overlap simulation are also 
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Figure 2: Power law and exponential fits for the mean multi-overlap tunneling time [r™"'i] 
in 3d and id. 



indicated in the figure. Both Pj{q) probability densities take their minimum at g = and 
we have pp^/pf^ > lo^ in 3d (L = 8) and pp^/pf'^ > lo^ in M {L = 8); compare also 
table 2. The improvement factors in computer time are directly proportional and close to 
these numbers which reflect the enhancements in visits of -Pj(O) as compared to canonical 
simulations. Multiplying the improvement factors with the average CPU times needed by 
the multi-overlap algorithm for a single realization on lattices of this size (see table |I|), it 
becomes clear that exploration of such barriers by means of a canonical MC simulation is 
simply impossible. 

We conclude this section with remarks about the performance of the algorithm and 
implications on the physics of the system. The multi-overlap algorithm flattens the free- 
energy barriers F^. If they were the only cause for the slowing down of the canonical 
dynamics, the multi-overlap autocorrelation time should be dominated by a random walk 
behavior between g = — 1 and q = +1 and scale proportional to A^ (in units of sweeps). 
Fitting the estimates of the mean autocorrelation time [r'^^'^jav, where the average is with 
respect to the realizations J', to the power-law form ln([r™"^]av) = a + z In(A^) gives 
z = 2.32 ± 0.07 in 3d and z = 1.94 ± 0.02 in Ad. The fits are depicted in figure 0. Their 
quality is bad, nevertheless they show that the slowing down is quite off from the theoretical 
optimum z = 1. Exponential fits ln([r™"^]av) = Cq + ciN are also depicted in the figure. 



Whereas in 3d the exponential fit is far worse than the power-law fit, it is the other way 
round in Ad. Hence, the smaller ^-value in Ad should not be taken seriously. 

The physically important conclusion is: the observed large autocorrelation times demon- 
strate that, in the model considered, canonical overlap barriers are not an exclusive cause 
for the slowing down of spin- glass dynamics below the freezing temperature. Therefore, r^ 
has to be a lower bound of the full canonical autocorrelation time r^'^^: 

r| < r^"" . (10) 

One should understand q as one relevant direction in a complex, multidimensional configu- 
ration space. By depicting free-energy barriers as function of q one projects on this direction 
and averages results over all other directions. 

3 Barrier Results 

We analyze our free-energy barrier densities relying on a variant of the cumulative distribu- 
tion function F. For a set of sorted data 

Xi < X2 < ■ ■ ■ < Xn (11) 



the (empirical) cumulative distribution function F{x), see for instance [^, is defined by 



i \ i 1 

-- — < F{x) <- + — for Xi<x< Xi+i , (12) 

n In n 2n 



where we use a straight-line interpolation in-between. Next we define a Q-tile |]3^ distribu- 
tion function as introduced in [^ 

p .,(F{x)ioiF{x)<0.5; 

^Ql-^J \ 1 _ ir(a;) for F{x) > 0.5 . ^ ^ 

This function peaks at the median Xmed of the data and takes there the value Fq = 0.5. For 
self- averaging data x the function Fq collapses in the infinite volume to 

0.5 for X = x ; 



Fn(x) , 

otherwise . 

Here x is the mean value. For non-averaging quantities the width of Fq stays finite. The 
concept carries over to observables which diverge in the infinite volume limit, when on each 
lattice size results are expressed in units of the respective median value, i.e. instead of an 
observable X the ratio x = X/X^^d is used. 
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Figure 3: Distribution function Fq (|T^) for the 'id overlap barriers (^ in units of their 
median value. 
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for the 4(i overlap barriers (^ in units of their 



3.1 Lack of self- averaging 

For the free-energy barriers (y) we have depicted our thus obtained Fq[Fq/ Fj^ ^^^ functions 
in figures |^ {?>d) and ^ {Ad). For each lattice the measured F% values were first sorted as 
function of J such that 

Fba < Fb,2 < ■■■ < F% , 
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Figure 5: Distribution function Fq ( ]TB| ) for the Sc? energies (^ in units of their median value. 
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Figure 6: Distribution function Fq ( ]I3| ) for the Ad energies (^ in units of their median value. 

where n is the number of realizations #^7 given in table |l|. Subsequently F| ^^^ was 
calculated as 



-^B,med 



^ [K,n/2 + ^B,i+n/2) (^ IS ^ven in our cases) 



and Fq computed for x = F'g/F^ ^^^. 

Both figures support that Fg is a non-self-averaging quantity. This is stronger in 4d than 
in 3d, because the inner lines belong in 3d to the larger lattices, whereas in Ad it is the other 
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Figure 7: Distribution function Fq (O) for the 3d barriers (ffl) of the multi-overlap algorithm 
in units of their median value. 



way round. However, in both cases there are marginal finite-size effects, whereas finite-size 
dependence of self-averaging is expected to be rather strong. This becomes obvious when 
comparing with an observable which is supposed to be self-averaging. Namely, figures ^ {3d) 
and ^ {4d) depict the same analysis for the internal energy (^. In 3d self- averaging of this 
quantity is obvious, whereas in Ad there is an irregularity when going from L = 6 to L = 8. 
As our simulation temperature in 4d is quite low, we think that this behavior is related to 
groundstate irregularities on small lattices (only the corresponding half of the distribution 
is affected). For both 3d and 4d the g-tile distribution function of the energy is strongly 
peaked around E / E^^d = 1, whereas the overlap barrier distributions are much broader. 

It is generally believed that, in contrast to the equilibrium autocorrelation times con- 
sidered here, non-equilibrium autocorrelations are self-averaging |l^. No sample-to-sample 
deviations have been reported for real experiments [0] and self-averaging is also used for 



measurements of non-equilibrium properties in MC simulations |Tl|, [T^. 

The multi-overlap algorithm eliminates the free-energy barriers which are visible in the 
Pj{q) probability densities. Let us therefore focus on the autocorrelations times of this 
algorithm and its barriers defined by 
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Figure 8: Distribution function Fq (O) for the Ad barriers (ffl) of the multi-overlap algorithm 
in units of their median value. 

We perform the analysis along our previous lines and show in figures {3d) and | {4d) the 
thus obtained Fq{F^^^ / F^^^^) functions. Lack of self- averaging is even more obvious than 
for Fg. In figure |^ {3d) there are (within the statistical accuracy) no finite-size effects visible 
and figure ^ {4d) exhibits a strong anti-self-averaging trend: Results for the larger lattices 
move to the outside instead to the inside. 



3.2 Finite-size scaling behavior 

In this final part of section ^ we discuss how data (experimental or MC) for non-self-averaging 
observables may be analyzed such that comparisons of results from different groups become 
possible. One has to investigate many samples and should report the finite-size scaling 
behavior for fixed values of the cumulative distribution function F (|12|). In particular 
this includes F = 1/2 which defines the median value. We exemplify this for the overlap 
autocorrelation time r^ (|^), but the method applies for non-self-averaging observables in 
general. 

From figures ^ and § it is obvious that the autocorrelation times r^ will have long tails 
towards large values. This implies that the mean value over all samples is a rather erratic 
quantity which is dominated by a few rare realizations. Table ^ collects the mean, median 
{F = 1/2) and maximum {F = 1 — l/(2r7,)) values for r^. The numbers in parenthesis indicate 
error bars in the last digits of the quantity given before. The results show that contributions 
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Figure 9: Fits ( ]TB| ) of the Sc? free-energy barriers Fg versus A^^/^ corresponding to the 
exponential finite-size scahng behavior ([16|) of r^. From down to up the hues are at 
16 F 



1,3,5,7,9,11,12,13,14, and 15. 



of the maximum values dominate to a large extent the mean values (just divide the maximum 
values by the number of realizations and compare the results with the mean values). The 
maximum values rely on realizations J of likelihood 1/n, what explains their very large 
errors. In contrast to the mean and the maximum, results for fixed F remain well defined 
as long as F stays away from its extreme values 1/n and 1 — \jn. In particular, note that 
samples with relaxation times too long to be measured can still contribute to determine the 
correct F values for smaller relaxation times. 

In the following we focus on our results for the free-energy barriers F^ at F = i/16 with 
i = 1, . . . , 15. For each F value we performed fits to the form 



F| = ai + as ArV3 



which corresponds to the exponential finite-size scaling behavior 



T-g _ pCii pa.2 A^^/^ 

r^ — e e 



(15) 



(16) 



suggested by investigations of autocorrelation times and barriers in the mean- field limit JIB, 
P^ |T^, |T^, ^. These fits are depicted in figures ^ and |TD|. Examples of the fit parameters Oi 
and 02 are collected in table ||; for all fits given there, the goodness-of-fit parameter Q |^ 
is smaller than 0.003. The average Q over all 15 fits is given in the figures. For consistent 
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Figure 10: Fits ( ]T3| ) of the Ad free-energy barriers F^ versus N'^^^ corresponding to the 
exponential finite-size scahng behavior ([16|) of r^. From down to up the hues are at 
16 F 



1,3,5,7,9,11,12, 13,14, and 15. 



fits the expectation for the Q-average is 1/2 and the quahty of our 3d and 4d exponential 
fits is unacceptable. We therefore try a power-law fit 



T 



B 



cN'^ 



(17) 



which corresponds to a fit of the form 



F'b 



ln(c) + a In(iV) 



These fits are depicted in figures |TT| and |1^. In 3d as well as in Ad the average Q-value is 
now almost perfect. Examples of the power-law fit parameters and Q- values are given in 
table 1^. They indicate that the distribution of Q- values is less perfect than their mean. Such 
uncertainties are an intrinsic limitation of MC simulations and become particularly severe 
when one is, as in our investigation, limited to rather small-sized systems. Having these 
limitations in mind, the over-all quality of the power-law fits is remarkably good. Our data 
favor them strongly over the exponential behavior. 

As function of F the exponent a = <y{F) varies smoothly and covers in Ad a range from 
0.8 at F = 1/15 to 1.3 at F = 15/16. In 3d the range is somewhat smaller, see table ^. Fits 
for F > 15/16 become erratic and it makes little sense to report them. A similar analysis 
for the autocorrelation times of the multi-overlap algorithm gives exponents a{F) which are 
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Figure 11: Fits (|TB|) of the 3d free-energy barriers F^ versus In(A^) corresponding to 
the power-law finite-size scahng behavior (|1^) of r^. From down to up the hues are at 
16 F = 1, 3, 5, 7, 9, 11, 12, 13, 14, and 15. 



5 


Average Q=0.50 




^^,„..^^^^ - 


4 




^,*''' 




3 




, 'S' ,- '" 




2 




_-""*" .--'"" 


--'-"" ,-"^ 




*'' --"'' „'"''' 


-'--*'' -'-"'' 


-""' ,--'""-'-'^ 


1 


'i'""""" ^ "-""-'-'- "'^-' 


:^-'5i;;;^^§^^^ 


''''''^^^ - 





^;;'-'-::;-'-;:--;HH:;;;;;;;:;;SH 


:::::f;:::t'::r;::::::::::::^= 


:'---::-'''"'I'Z^'^ 




• ;:..-;;:;;;:;;-;;:;;;::;;;;;;;;;:::;;: 


-:::'--::::$--'■';;;,-'—-''" 




9 

o 


IliigS--"""' 




" 



5.5 



6.5 



7 
ln(N) 



7.5 



8.5 



Figure 12: Fits (|TBp of the 4(i free-energy barriers Fg versus In(A^) corresponding to 
the power- law finite-size scaling behavior (|1^) of r^. From down to up the lines are at 
16 F 



1,3,5,7,9,11,12,13,14, and 15. 



larger. 



a^^\F) ^ a%{F) + 1 . 



This re- iterates and sharpens our previous observation that relevant barriers exist, which are 
invisible in the overlap variable q. 
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4 Summary and Conclusions 

We have investigated free-energy barriers in the Parisi order parameter (|T|). The results are 
sample dependent and non-self-averaging on the (admittedly rather small) simulated systems. 
The power-law behavior ([TtD of the Markov autocorrelation times r^ as defined in eq. (H) is 
favored over the exponential behavior (0). To the extent that this behavior extrapolates to 
the infinite volume limit and that our methods relate to those of Ref. [|l^, |l^, |18|, [19|, ^, it 
means that both 3d and 4d are quite far away from the d —>■ oo mean-field theory limit. As 
relevant barriers are still found in the autocorrelations of the multi-overlap algorithm, such 
a relation is far from clear. 
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Tables 



Table 1: Statistics: Number of realizations #J', average number of mega-sweeps per 
realization ngw O'lT-d average single 375 MHz processor CPU time per realization in hours 
(h) or seconds (s) as benchmarked on the CEA T3E. 





3d 


4d 


L 


i^J 


'^sw 


CPU 


i^J 


^sw 


CPU 


4 


8192 


0.2 M 


6.32 s 


4096 


0.4 M 


76.6 s 


6 


8192 


1.0 M 


113 s 


4096 


3.7 M 


1.02 h 


8 


8192 


7.6 M 


0.54 h 


1024 


49.3 M 


42.66 h 


12 


640 


154.0 M 


36.97 h 









Table 2: Mean, median and maximum values for the Markov autocorrelations times Tg ([^. 



L 


4 


6 


8 


12 


3d 


Mean 


61(29) xlO^ 


10(06) xl06 


56(45) xlO^ 


13(10) xlO^ 


3d 


Median 


237(05) xlO^ 


690(02) xlO^ 


152(04) xl02 


444(05) xl02 


3d 


Maximum 


22(20) xlO^ 


32(06) xlO^ 


35(33) xlOi° 


53(30) xlO^ 


Ad 


Mean 


94(34) xlO^ 


23(08) xlO^ 


26(23) xlO^ 




4d 


Median 


807(18) xlO^ 


379(11) xl02 


117(07) xlO^ 




Ad 


Maximum 


13(11) xlO^ 


21(07) xlO« 


22(21) xlO^ 





Table 3: Fit parameters for the free-energy barriers fits j^TSj); Q < 0.003 for all cases. 





3d 


M 


F 


ai 


0-2 


a-i 


(22 


1/16 


4.59(2) 


0.424(03) 


6.30(02) 


0.274(03) 


4/16 


5.13(3) 


0.442(04) 


6.77(03) 


0.276(03) 


8/16 


6.38(4) 


0.382(05) 


7.13(05) 


0.301(06) 


12/16 


7.39(6) 


0.476(11) 


7.77(11) 


0.344(10) 


15/16 


9.71(9) 


0.538(14) 


8.85(12) 


0.435(15) 



Table 4: Free energy barriers fits j^TSj): Fit parameters and goodness 



of fit. 





3d 


Ad 


F 


He) 


a 


Q 


ln(c) 


a 


Q 


1/16 


2.80(03) 


0.830(06) 


0.04 


3.58(05) 


0.804(08) 


0.05 


4/16 


3.30(04) 


0.857(08) 


0.77 


3.68(06) 


0.860(08) 


0.13 


8/16 


4.10(06) 


0.883(10) 


0.52 


3.68(11) 


0.958(18) 


0.71 


12/16 


5.37(11) 


0.930(20) 


0.51 


3.71(22) 


1.105(32) 


0.81 


15/16 


7.35(14) 


1.075(28) 


0.02 


4.37(26) 


1.302(42) 


0.94 
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